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Abstract 

We consider differences between logr(x) and truncations of certain clas¬ 
sical asymptotic expansions in inverse powers of x — A whose coefficients are 
expressed in terms of Bernoulli polynomials B n { A), and we obtain conditions 
under which these differences are strictly completely monotonic. In the sym¬ 
metric cases A = 0 and A = 1/2, we recover results of Sonin, Norland and 
Alzer. Also we show how to derive these asymptotic expansions using the 
functional equation of the logarithmic derivative of the Euler gamma function, 
the representation of 1/x as a difference F(x + 1) — F(x), and a backward 
induction. 


1 Introduction 

The Euler gamma function T(x) is widely regarded as the most important special 
function. Good accounts and many formulas for T(x) can be found in [4], [5, §6.3], 
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[6, Chapter 1], and [12, Chapter 12]. Various bounds are known for T(x) and if(x), 
the logarithmic derivative of T(x), also called the digamma function, e.g. [1], [2], 
[5, 6.3.21], [7], [9, Appendix C], [10]. The logarithm of the gamma function has the 
classical asymptotic expansion 

(1) log T(x) ~ log V2tt + (x ~ |) log(x — A) — (x — A) 

+ - a) 1 -, 

^ n(n- 1) 

which we consider for real x > A and A G [0,1]. As usual, B n ( A) denotes the nth 
Bernoulli function. Our focus here is on differences (and their derivatives) between 
logT(x) and the truncations 

L n ( A; x) := log V2tt + (x — §) log(x — A) — (x — A) 

n =2 v ' 

of this asymptotic series. In the most basic case, we are interested in conditions under 
which these truncations provide strict upper or lower bounds. These inequalities 
together with the functional equations provide an effective method for calculating 
logT(x) and its derivatives with arbitrary accuracy. 

Example 1.1. A pleasant property of (1) is that we can differentiate both sides to 
obtain asymptotic expansions for the derivatives of logT(x), that is, the derivatives 
of the truncations A;x) provide asymptotic expansions for the derivatives of 
logT(x). For instance, in the particularly symmetric case A = 1/2, we have the 
asymptotic expansion 

(2) i,(x) ~ log(* - 1/2) - f; g2n / /2) (x - 1/2)- 2 3 4 ”. 

n= 1 

Our main theorem implies the following result of Norlund [11, §56, p. 106]: if the 
series in (2) is truncated at n = N with N > 0 and even, the resulting sum is a 
lower bound for if(x) valid for all x > 1/2; if N is odd, it is an upper bound. In 
other words, successive truncations of (2) provide lower and upper bounds valid for 
all x > 1/2. The first examples of these bounds are 

(3) il](x) > \og(x — 1/2) 
and 

(4) ip(x) < log(x — 1/2) + (x — l/2) _2 /24 . 
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A function / is said to be completely monotonic on an interval or half-line / 
if, for all x G I and all integers n > 0, the derivative f^ n \x) exists and satisfies 
(—l) n /( n )(x) > 0. If this inequality is strict, that is, if (—1 ) n f^ n \x) > 0, then 
/ is called strictly completely monotonic on I. For an introduction to completely 
monotonic functions and their properties we refer to [13, Chapter IV], 

Our main result, which is proved in Section 4, is the following. As indicated in 
Example 1.4, this is an extension of a result of Alzer [1, Theorem 8], which concerns 
the case A = 0. 

Theorem 1.2. Let N > 1. Suppose that either N is even and A G [0, Ao] where Ao 
is the unique root of A) in [0,1/2], or N is odd and A G [Ao, 1/2] where Ao is the 
unique root of Bn+i(\) in [0,1/2], Then, 

(_l) rAr / 2 l (log r(x) - L n (\-x)) 

is strictly completely monotonic on (A,oo). 

We note that Ao < 1/4 and that Ao —* 1/4 very rapidly as N increases. This 
is explained in Example 3.3, which also illustrates that Theorem 1.2 would be false 
without the conditions on N and A. On the other hand, the restriction to the case 
A G [0,1/2] is made for the sake of exposition, and we invite the interested reader to 
adjust the details for the case A G [1/2,1], 

Example 1.3. The particular case A = 1/2 and N = 3 shows that 

logT(x) - log Vrii - [x — \) (log [x - |) — l) + Jkjd_ 

x — 2 

is strictly completely monotonic on (1/2, oo). This generalizes the inequality (4). 

Example 1.4. In the case A = 0, we find that the functions if(x) — L2m(0;x), 
explicitly given by 

M 

if(x) - log - (x - |) log(x) +X-J2 2n(2n-lW^ 

are strictly completely monotonic on (0, oo) if M is even. Likewise, the functions 
L 2 m(0; x) — ip(x) are strictly completely monotonic on (0, oo) if M is odd. This is a 
result of Alzer [1, Theorem 8]. 

This paper is organized as follows. In Section 2, we give a first illustration of 
our approach by proving the claims made in Example 1.1. These bounds for the 
digamma function are then generalized in Section 3, before we prove the main result 
in Section 4. A very brief indication of how to use these bounds for numerical 
calculations is then given in Section 5. 
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2 Proving bounds for the digamma function 

In this section, we illustrate our approach by proving the claims made in Example 1.1. 
In particular, we show, by a method we believe new, how to derive asymptotic 
expansions for such as (2), whose successive truncations provide lower and 

upper bounds. 

Assuming there is an approximation of the form ^(x) ~ F(x), then the functional 
equation -0(x + 1) = ijj{x) + 1/x would give 

(5) F{x + 1) — F{x) « 1/x. 

Our method is to find functions F giving upper (resp. lower) bounds in (5) and then 
show by a “backward induction” that these functions provide bounds for i^(x). We 
use one additional fact about the digamma function iJj(x), namely the asymptotic 
formula 

(6) ^(x) = logx + o(l), x —* oo. 

For completeness, we prove this at the end of the article. 

Example 2.1. To illustrate our method, we first obtain the lower bound (3). By 
convexity, 

f x+l / 2 d u 1 

log(x + 1/2) - log(x - 1/2) = / — >-. 

Jx- 1/2 U X 

This inequality and the functional equation give 

'ijj(x) > ij){x + 1) — (log(x + 1/2) — log(x — 1/2)) 
for all x > 1/2, or, equivalently, 

(7) ^(x) — log(x — 1/2) > ij){x + 1) — log(x + 1/2). 

Now we prove (3) by an induction: repeating (7) with x replaced by x + 1, x + 
2, ... and using (6), we see that 

ip{x) — log(x — 1/2) > liminf — log (y — 1/2)) = 0. 

y—>-oo 

The key to our treatment of the general case is to develop a suitable expression 
for 1/x with the aid of the Euler-Maclaurin sum formula. The following version of 
Euler-Maclaurin, given in [9, Appendix B] is convenient for our application. 
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Theorem 2.2 (Euler-Maclaurin). Let K be a positive integer, and suppose that f 
is K times continuously differentiable on [ a,b ]. Let B k f) denote the k-th Bernoulli 
function. Then, 


f( n ) = / f( x ) dx 

a<n<b J a 

+E - s t ({«})/ < ‘- 1 ) (a)) 

k =1 

B K({x})f (K \x)dx, 

where {t} := t — |_£j. 

We now show the following result of Norlund [11, §56, p. 106] stated in Exam¬ 
ple 1.1. 


Theorem 2.3. The asymptotic expansion (2) for if (x) holds as x —> oo. Moreover, 
if the series in (2) is truncated at n = N with N > 0 and even, the resulting sum is 
a lower bound for if(x) valid for all x > 1/2; if N is odd, it is an upper bound. 


Proof. Take f{t) := l/{x+t), a = —1/2, b = 1/2 in Theorem 2.2. Because B k { 1/2) = 
0 for all odd k, we restrict the series to even indices. Starting with the one-term (!) 
sum 


1 

x 


E 


—l/2<n<l/2 


X + Tl 


we have 


1 

x 


r 1/2 d t 

1-1/2 X + t 


\ ^ B 2n (l/2) ( 1 _ 1 \ td 

^ 2 n v (x + l/2) 2n (x-l/2) 2 "7 N 


where 

f 1/2 d t 

For the claimed inequalities, we establish the sign of Rn with the aid of several 
properties of Bernoulli functions, all of which are given in [9]. For —1/2 < t < 0, 


B 2 N({t}) — B 2 n( 1 + t)— B 2 N(—t). 
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Thus, we have 


Rn — ~ 


1/2 

B2N(t) 



1 _ 1 \ 

t) 2N+1 (x + t) 2N+1 ) 


d£. 


Next, 


B2N+i(y) — (2iV + 1) / B 2 N{t) d t. 


Integrating by parts and noting that the integrated term vanishes because T> 2 jv+i(1/2) 
B 2 n+ i( 0) = 0, we find 


,1/2 


Rn — 


B 


2N+1 


H) 


.(; X-t ) 2N + 2 ( x + t)27V+2 

For 0 < t < 1/2, sgn(ij 2 V+i (f)) = (—1) 7V+1 . Also, in this range, 

(x-t)- 2N ~ 2 - (x + t)~ 2N ~ 2 >0 


d t. 


by monotonicity. Thus Rn < 0 for N even and Rn > 0 for N odd. 

Setting 

(8) F n (x) := log(x - 1/2) - £ B2 "£/ 2) ~ 1 / 2 

n =1 

we find 

-0(x + 1) - ^(x) = - = Fjv(x + 1) - Tjv(V) + R N - 

X 

Thus 

-0(x) - F n {x) > Tp(x + 1) - F N {x + 1) 
for N even (resp. < 0 for N odd). 

Assuming that N is even, we proceed as in Example 2.1. We find inductively 
that, for all x > 1/2, 

i/>{x) - F n (x) > liminf(^(j/) - F N (y)) = 0, 

y—> oo 

which proves the lower bound for ip{x). An analogous argument shows that i(j{x) — 
F n (x) < 0 for N odd. □ 
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Example 2.4. Integrating formula (2) termwise, we find, for a suitable constant C, 


logT(x) ~ C+(x — 1/2) log(x — 1/2) — (x — 1/2) 


^2n(l/2) 

2n(2n - 1) 


l/2) l ~ 2n . 


We next exponentiate and then determine C by the well-known device of using the 
duplication formula for the gamma function to connect T( 2x) with r(x + 1/2) and 
T(x). Letting x —> oo, the terms in the sum go to zero, and we find C = log \[2rx. 

As a consequence of Theorem 2.3, we obtain that, for any positive integer N and 
all x > 1/2, 


(9) 


r 2 w-i(z) < r(x) < t 2 n(x), 


where 


T N (x) = V2n (x - l/2) x ~ 1/2 e~( x ~ l/2) 


x exp l ^ 

^ 71 = 1 


B 2n (l/2) 
2n{2n — 1) 


(x - 1/2) 


1-2 n 


These bounds for the gamma function are given by Norland [11, §58, p. Ill], who 
attributes them to Sonin. 

Note that the inequalities in (9) are opposite from those of Theorem 2.3, because 
the integration has changed the signs in the summations. 


3 Generalized bounds for the digamma function 


Let A G [0, 1]. As x —> oo, we have also the asymptotic expansion 
(10) V-W ~ log(x - A) - f)(-l)”AW (s - A)-” 

71 

n =1 


This generalization of the formula (2) in Theorem 2.3 can be derived by our method. 
It appears also, e.g., in [8, §2.11, p. 33]. Note that here, except at A = 0,1/2,1, the 
Bernoulli functions with odd indices are nonzero. For A = 0, B n ( 0) = B n: the nth 
Bernoulli number, and (10) becomes the well-known asymptotic expansion [5, 6.3.18] 

1 °° D 

(n) ~ i°g x - — - Y x ~ 2n ■ 

n= 1 
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At first glance, this formula looks very similar to the expansion (2). However, in 
addition to having an extra term and a different center of expansion, the coefficients 
of the last series are Bernoulli numbers rather than Bernoulli functions evaluated at 
1/2; the connection between these is given by the formula B n ( 1/2) = —(1 — 2 l ~ n )B n 
for n = 0,1,... [5, 23.1.21], 

Denote the truncation of the asymptotic series in (10) at n = N by 



( 12 ) 


n= 1 


The following generalization of Theorem 2.3 specifies conditions under which these 


approximations to 0(x) provide lower or upper bounds that are valid for all x > X. 
As in the case of our main result, Theorem 1.2, which is a natural extension of these 
bounds, the restriction to the case A G [0,1/2] is made for expository reasons. 

Theorem 3.1. Let N > 1. Let Ao be the unique root of B^(X) in [0,1/2] if N is 
even, and the unique root of B^+ i(A) in [0,1/2] if N is odd. 

• Suppose that N = 1 modulo 4 and A G [A 0 ,1/2]. Then, we have 0(x) > 
Fj\t(X] x) for all x > X. If N = 3 modulo 4, then the inequality for 0 is 
reversed, that is, if(x) < F^(X]x). 

• Suppose that N = 2 modulo 4 and X G [0, Ao]. Then, we have 0(0) > Fn( A;x) 
for all x > X. If N = 0 modulo 4, then the inequality for 0 is reversed. 

We note that Ao < 1/4 and that Ao —> 1/4 very rapidly as N increases; see 
Example 3.3. 

In preparation for the proof of Theorem 3.1, we deduce the following approximate 
version of the functional equation 0(x + 1) — 0(x) = 1/x. 

Proposition 3.2. Let N > 1. Then 



Proof. As in the case A = 1/2, we set f(t) := 1 /(x+t) and apply the Euler-Maclaurin 
summation formula to the one-term sum 
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We find 


1 

x 


r x j- BJ1 ~ x l 

x + t ^^n\(x+l — X) n 



s N m 

(x + t)^ 1 


Bn{ 1-A) \ 
(x - A) n ) 


and it remains only to note that B n (l — A) = (— l) n B n (X). □ 

Proof of Theorem 3.1. First, assume that N = 1 modulo 4 and that A G [Ao, 1/2]. In 
light of Proposition 3.2 and the backward induction already used in the case A = 1/2, 
it suffices to show that 


U-A 


Rn '■= 


M W) 

I- A (x + t) N+l 

To determine the sign, we split this integral as 

p _ p/ | P " ._ [ X B N({t}) i, 

R n -R+R .-J^ {x + t)N+1 dt + 


dt > 0. 


“1—A 


_MW) 

(x + t) 


N +1 


dt. 


For t e (—A, 0), we have -Bjv({t}) = B N (l + t) = —B N (—t ) since N is odd, and hence 


R' = 


B N (-t ) 

/_A (x + t ) iV+1 


-Bjv(t) 


dt + 


-Bjy(t) 

(x + t)^ 1 

1 


dt 


(z -t) 


JV+1 


(x + t) 


JV+l 


dt. 


Integrating by parts and observing that the boundary terms vanish, we find 


R'= / (B N+1 (t) - B n+1 ( A)) 


+ 


(x — t) Ar+2 (x + t) N+2 


dt. 


Likewise, splitting R" according to J A ' A = f^ 2 + /p 2 A , then integrating by parts, 
we obtain the counterpart 

R " = l {BN+l{t) ~ B ' v+l(A)) ((i + i)«« + (* + 1 - i) K+2 

In order to bound R^ — R' + i?", we determine, in both integrals, the signs of each 
of the two factors. Since N = 1 modulo 4, the polynomial B N+ i{t) is decreasing for 
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0 < t < 1/2 [9, (B.12)]. It follows that B^+iit) — Bn+i(X) > 0 for 0 < t < A, and 
B]y + i(t) — B n+1 (X) < 0 for A < t < 1/2. Since the function 

1 1 

(x — t) N+2 (x + t) N+2 

is increasing and positive for t such that 0 < t < x, we find that 

2 


R' > 


N+2 


X 


{B n+ i(t) — B N+ i(X))dt. 


On the other hand, for all t G [0, 1], 

1 


< 


(x + t) N+2 (x + 1 — t) N+2 X N+2 ’ 
which allows us to bound 


' 1/2 


R" > 


N+2 


X 


( B]y + i(t ) — BN + i(X))dt. 


Combining the two inequalities, we obtain 

2 


Rn > 


x 


N+2 


1/2 (B N+1 (t) - B N+1 (X))dt = Bn +^ X) 


x 


N+2 


Since A 0 is the zero of B N+ i(A), which is decreasing on [0,1/2], and since A > A 0 , 
we conclude that Rjy > 0. It is straightforward to adjust this proof to cover also the 
case N = 3 modulo 4. 

Next, assume that N = 2 modulo 4 and that A G [0, A 0 ] • Again, it suffices to 
show that 


Rn 


-dt > 0. 


I- A i x + t ) N+1 

Let £{t) be the linear function whose graph intersects the graph of l/(x + t) N+1 at 
t = Ao and t — 1 — Ao- Since 5 at({/}) < 0 for Ao < t < 1 — Ao, and > 0 for 

—A 0 < t < A 0 or 1 — A 0 < t < 1 + A 0 , it follows that 

1 




-m 


> 0 


_(x + t) Ar+1 

for all t G [—A 0 ,1 + A 0 ], with strict inequality if t ^ A 0 and t ^ 1 — A 0 . In particular, 


writing £(t) = mt + b, with m < 0, we have 

/>1—A 


-1-A 


(13) 


Rn > 


I- A 


BN({t})(mt + b)dt = m / tB N ({t})dt. 


’ —A 
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The last equality is a consequence of the fact that the average value of Bjsr({t}) over 
a period is zero. Splitting the last integral as f^ x ' X = f°_ x + A , we obtain 



tB N ({t})dt 


fl-A 


tB N ( 1 + t)dt + 


hEbv(t)df. 


'-A 


After a change of variables in the first integral and the observation that tB^(t)dt = 
0 , this simplifies to 


'1—A 


'-A 


hEbv({f})dt = / tBN(t)dt — / Bjy(t)dt — — I 5jv(t)dt < 0, 


' 1 —A 


1—A 


Since m < 0, we conclude from (13) that Rn > 0. Again, it is straightforward to 
adjust the proof to cover the case when N = 0 modulo 4. □ 


Example 3.3. The equality of the generating functions 

B / ] x x 2n 1 / xe x / 4 xe~ x ! A A 1 xe 3 ^ 4 yA (|) x n 

2n viJ (2n)! ~~ 2 Ve* - 1 ~~ e~ x - 1 J ~~ 2e x / 2 - 1 ~~ ^ 2 n 7d 

implies that £> 2 at(1/4) = 2~ 2A, T> 2 jv(l/2). Since sgnS 2 Jv(l/4) = sgn B 2j v(1/2) and 
£> 2 at(A) is monotonic on (0, 1/2), we conclude that the unique root Ao of the Bernoulli 
polynomial B 2 n( A) in [0, 1/2] satisfies Ao < 1/4. Theorem 3.1 therefore implies that, 
for x > 1/4, 

V>0) > T\(l/4; x) = log(x - \) - -^4 , 

x- ^ 

V>(x) < F 3 (l/4 ;x) = log(x - A) - -^4 + ^ » 

x - 4 (a; - iY (x- j) 6 

and so on. On the other hand, for instance, the difference ip{x) — F 2 (l/4; x) changes 
sign at some x in the interval (1/4, oo). 


Example 3.4. Applying Theorem 3.1 in the special case A = 0, we find that trun¬ 
cating the classical asymptotic series ( 11 ) results in strict upper and lower bounds 
for il>(x), valid for all x > 0. The first two such bounds, namely 


v-w > logW - T - J-j 


- 0 (x) < log(x) 


1 

2 x 


1 


+ 


12 x 2 120 x 4 

are proved also in [7, Theorem 5]. That these inequalities hold more generally is 
implied by [1, Theorem 8 ]. 
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4 Complete monotonicity 


We are now ready to generalize onr approach to give a proof of Theorem 1.2, our main 
result. First, we note that, combining previous observations, it is straightforward to 
make the error term in the asymptotic expansion (1) explicit. The special case A = 0 
of (14) appears, for instance, in [3, Theorem D.3.2], and equation (15) is given in 
[11, §53, p. 101, (6)]. 

Proposition 4.1. Let N > 1. Then for x > X 

(14) logr(l) " iA ' (A;i) = _ v/”pTp di - 

Proof. This follows by integration, and the fact that logT(x) — L N { A; x) converges 
to 0 as x —> oo, if we can show that 

( 15 ) ^( x ) ~F n (X;x) = J ^ dt ’ 

where F/v(A; x) is dehned as in (12). Note that Proposition 3.2 can be restated as 

fj(x) - F n ( A; x) = if{x + 1) - F n (X] x + 1) + J dl 

Iterating this relation, we hnd that, for all integers n > 1, 

if(x) - F n (X] x) = if{x + n) - F n (X] x + n) + [ dt, 

J -A \X + t) 


and equation (15) follows from taking the limit n —> oo. 


□ 


Proof of Theorem 1.2. Consider, for illustration, the case that N = 1 modulo 4. 
Under the given assumptions, Theorem 3.1 states that f(x ) = Ljy( A;x) — logT(x) 
satisfies f'(x) < 0. Integrating as in Example 2.4 further shows that /(x) > 0. 

We now indicate why, more generally, (—l) n / (n )(x) > 0. Note that (14) implies 


that 


/ 0 ) - /0 + 1 ) 


1 f^Bsdt}) 
Nj_ x (x + t)» m 


Differentiating this relation, we observe that 


/<”>(!)- /<")(!+ 1 ) 


iy ,w. 

J N 



Mid* 

(x + t) N + n ’ 
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where (N) n = N(N + 1) • • • (N + n — 1) >0. By the usual backward induction, it 
follows that (—1 ) n f^ n \x) > 0 , if we can show that 



£*({*}) 

(x + t) N+n 


dt > 0 . 


The proof of Theorem 3.1, which shows this inequality in the case n — 1, extends in 
a straightforward manner to show that these inequalities hold in general. □ 


5 Numerical calculations 


It is generally infeasible to make a calculation of prescribed accuracy using an asymp¬ 
totic expansion. We can do this for ip(x) (as well as, more generally, for logT(x) and 
all of its derivatives), however, because of the following two special properties. 

First, we have explicit error terms for the series (2). The corresponding extensions 
to logT(x) and its derivatives (compare Theorem 1.2) as well as to values of A other 
than A = 1/2 are straightforward. 

Corollary 5.1. Let F N (x ) be defined by ( 8 ) and suppose x > 1/2. Then \ip(x) — 
F/v(x)| is smaller than the first omitted term of the series, and, using the average of 
consecutive terms, 


fj(x) 


F n (x) + F n+ i(x) 
2 


^ 1 | L$2N +2(1/2) | 1 /r) ^ _2JV—2 

S 2 21V + 2 /j 


Proof. ip(x) lies between F^(x) and Fn+ i(x). Thus its distance from F^(x) is smaller 
than \F n+1 (x) — Fjy(x) |, the absolute value of the last term of Fjv+i(x). Further, 
the distance of if(x) from the average of the two partial sums is at most half their 
separation. □ 


Second, each fixed partial sum of an asymptotic series of the type ( 2 ) is more 
accurate at larger arguments, and by applying the functional equation repeatedly, 
we can approximate ip at a larger argument. 

As an illustration, we find ^(1) (which equals — 7 ) to 20 decimal places of accu¬ 
racy. For a number I\ of our choice, write 

V>( 1) = tp{K) - 1 - 1/2 - 1/3- 1 /{K - 1 ), 

and then approximate ip{K). Say we take K = 100 and use the first four terms of 
the series in (2). We find 

V 24 7 / 960 31/8064 127/30720 

° g ' + 99732 - 99.54 + 99.56 ^58 
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as an approximation to 0(100). The error is guaranteed to be smaller than the first 
omitted term of the series, 


(511/67584) 99.5 -10 < 10 " 22 


We find 0(1) « —0.57721566490153286061. (The analogous calculation using the 
first four terms of the series in (11) yields the same result.) On the other hand, the 
most accurate approximation to 0(1), using the fixed choice K = 100, is achieved 
when we use, roughly, the first 300 terms of the asymptotic series. In that case, the 
error is less than 10 -2 ' 2 . 

6 The asymptotic formula ( 6 ) 

If we logarithmically differentiate the Weierstrass product formula for T(x) [5, 6.1.3] 
(unfortunately listed there under Euler’s name), we get 



From this we see that if M is a positive integer and M — x = 0(1), then 0(x) = 
0(M) + 0(1/M). To prove ( 6 ), it suffices to show 


(16) 


4>(M + 1) = log A/ + 0(1/M). 


Let L be a large positive integer and write 



and group terms in blocks of length M. We find upon rearranging that 


Slm = Y-- y — = log M + 7 + 0(1/M) + 0(1/L). 
f n n 


m 1 (L+1)M 


n= 1 n=LM+l 


Letting L —> oo, (16) and hence ( 6 ) is proved. 
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